Spin orbit interaction and zitterbewegung in symmetric wells 



Esmerindo Bernardes,'[] John Schliemann, 2,3 J. Carlos Egues, 1 ' 4 ' 5 ^ and Daniel Loss 3 ' 5 

'institute/ de Fisica de Sao Carlos 
Universidade de Sao Paulo 
Av. do Trabalhador Sao-carlense, 400 CP 369 
13560.970 Sao Carlos, SP, Brazil 
2 Institute for Theoretical Physics 

University of Regensburg 
D-93040 Regensburg, Germany 
3 Department of Physics and Astronomy 
University of Basel 
CH-4056 Basel, Switzerland 
4 Department of Physics and Astronomy 
University of Basel 
CH-4056 Basel, SwitzerlandK 
5 Kavli Institute for Theoretical Physics 
University of California, Santa Barbara, 93106 California, USA 
(Dated: February 1, 2008) 

Recently, we hav e introduced a nov el inter-subband-induced spin-orbit (s-o) coupling [Phys. Rev. Lett. 



99, 076603 (2007); :ond-mat/0607218] arising in symmetric wells with at least two subbands. This new s-o 
coupling gives rise to an usual zitterbewegung - i.e. the semiconductor analog to the relativistic trembling 
motion of electrons - with cycloidal motion without magnetic fields. Here we complement these findings by 
explicitly deriving expressions for the corresponding zitterbewegung in spin space. 



1. INTRODUCTION 

A controllable coupling between the spin and the orbital degrees of freedom in nanostructures is a highly desirable ingredient 
in the emerging fields of semiconductor spintronics and spin-based quantum computation and communication [JI]]. The gate- 
tunable [[]] Rashba s-o interaction [[J present in two-dimensional electron gases (2DEGs) formed in structurally asymmetric 
confining potentials offers such a possibility. The seminal proposal of Datta and Das [Q] of a spin field effect transistor highlights 
the use of the Rashba spin-orbit interaction to coherently rotate spins. More recently, the Rashba s-o interaction has been 
proposed as a convenient coherent mechanism for manipulating entangled electrons ("flying qubits") in the solid state [ ]. 
Interestingly, systems with bulk inversion asymmetry exhibit an additional spin orbit interaction, termed the Dresselhaus s-o 
whose interplay with the Rashba s-o can be used to devise a robust spin field effect transistor [[]]. 

So far, s-o effects in 2DEGs have only been investigated in (asymmetric) quantum-well systems with a single confined state. 
Recently, we have introduced a new type of spin-orbit interaction present in quantum wells with two confined subbands [[]]. 
Unlike the ordinary Rashba s-o interaction, ours is non-zero even in symmetric wells - as it arises from the inter-subband- 
coupling between the lowest (even) and the first (odd) excited states of the well. As pointed out in Ref. [Q, this new inter- 
subband-induced s-o coupling leads to an unusual dynamics of injected spin-polarized wave packets - the zitterbewegung or 
trembling motion - with cycloidal trajectories without magnetic fields. This is qualitatively different from the zitterbewegung 
in the presence of Rashba or Dresselhaus s-o interactions [[]]. Here we complement the investigation in Ref. [Q by deriving 
explicit expressions for the zittebewegung in spin space, i.e., the spin dynamics of the injected electron. In what follows, we first 
present our new Hamiltonian with its eigenvalues and the corresponding time-evolution operator, and then proceed to determine 
the spin dynamics of injected electrons. 

2. EFFECTIVE HAMILTONIAN: NOVEL SPIN-ORBIT TERM 

Starting from the usual 8x8 Kane model one can derive a 2 x 2 effective Schrodinger equation for the conduction electrons by 
a folding down process. We have performed this procedure [Q] for a symmetric quantum well with two subbands with edges at 
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the quantized energies s e (lowest level) and s„ (first level). In terms of the real-spin cr and pseudo-spin r (describing the subband 
degree of freedom) Pauli matrices, our 4x4 Hamiltonian for an electrons with effective mass m* can be cast in the compact form 



En ± e„ 



\ 2m* / n v ' 2 



where the new inter-subband-induced s-o coupling r\ is given by 



V = - 



F 2 



(e s + a)- 



p2 

— (e\d z V(z)\o) + 



6y 
El 



(E g + A) 



^{e\d z h(z)\o), 



(1) 



(2) 



where the Kane matrix element is P — -iH(S \p x \X)/mo, with \S ) and \X) being the band-edge (F point) periodic Bloch functions, 
the fundamental and split-off band gaps in the well region are E g and A, respectively, \e) and \o) denote the confined ground 
and first excited states of the well, V(z) is the Hartree-type contribution to the potential, h(z) is the profile function defining the 
quantum well, and 8 V , <5a are the offsets between the valence bands in the well and barrier regions. 

Note that "H bears a close similarity to the usual Rashba Hamiltonian [Q; however, the corresponding s-o coupling r\ can be 
non-zero even for symmetric wells - this can be easily understood by noting that r\ contains matrix elements between the ground 
and excited states. For a symmetric square-well potential, for instance, the first term in Eq. (||) vanishes (i.e. there is no Hartree 
contribution) and the second term reduces to 



2P 2 



T] = 



E 2 



(E g + A) J 



<p (a)(p e (a), 



(3) 



where <pi(z) i — e,o denote the confined well wave functions. Note that d z h(z) = -6{z - a) + 6(z- a) since h(z) = ®(a— z) + ®(z— a) 
[0(z): Heaviside function] defines a square well of width L = 2a. As shown in Ref. ([[]]) the magnitude of r\ can be comparable 
to that of the usual Rashba s-o coupling. 

The Hamiltonian (p can be easily diagonalized yielding the eigenvalues 



fi ± = Ek + SO, 



(4) 



where = ir- f + e + and ftQ = ^e 2 + rfk^. 



3. TIME EVOLUTION OPERATOR 



The time evolution operator U = e\p(-i ( Ht/K) can be straightforwardly obtained from Eq. (jl|). In the plane wave basis 
{|fc M , cr-.)}, we find 



U = e- iet,/h jcos(Qf)l 91 + i [e_r : <g> 1 - j]T x ® (k x <r y - k y cr x )] S1 ^ f) 
Next we use U to calculate the time evolution of the spin operators. 



(5) 



4. ZITTERBEWEGUNG IN SPIN SPACE 



To determine the spin dynamics in our system we have to calculate the time evolution of the tensor product 1 ® cr 1 , i = x, y, z. 
In the Heisenberg picture we have cr' H (t) = t/ 1 T ® cr'U, where cr"s denote the Pauli matrices at t = 0. For the z component we 
find 



(6) 



ai,(t) = { cos(Qf)l ® 1 - i \e-T ® 1 - tjt x ® (k x cr y - kyCT*)] 

{ L 1 nil 

l®cr |cos(Q01 ® 1 + i [e-r ® 1 - 7]T X ® (k x cr y - k y cr x )j & 



Tl 2 

cr^(f) = 1 ® a 2 - - -t-r x ® (k x a- x + k y cr x ) sin(2Qf) - - [e-T]T y ® (k x cr x + k y cr x ) + J] 2 k 2 j sin 2 (Of) 



n 2 n 2 



(7) 
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Similarly, for the x and y components we obtain 

cri(f) = 1 ® cr 1 + ^r r <g> sin(2Q0 + [2e_r v ® cr ; - 77I ® (Jt^tr* + k y cr x )\ sin 2 (Qf), (8) 
nil n iir L J 



cr^(f) = 1 <g> cr- v + ^-V ® cr' sin(2Qf) + [2e-r v <g> cr - 77I ® (fc x cr r + ^,cr x )] sin 2 (Qf). (9) 

In deriving the above equations we have used the algebraic relations (cr 1 ) 2 - \, i - x,y,z and cr'cr-' = -cr'cr' = icr* 7, A: = x, y, z 
(or any cyclic permutation), obeyed by the Pauli matrices. The r\i - x, y, z, also obey these relations. Equations describe 
the spin dynamics of injected electrons due to the s-o coupling 77. Together with the equations for xn(t) and ynit) derived in 
Ref. [[J, we have now the complete dynamics of an injected electron - i.e., the time dependence of both the spin and orbital 
degrees of freedom. As we discuss next, this dynamics are intrinsically linked due to the spin-orbit coupling and the fact that the 
z component of the total angular momentum j is conserved. 



5. SPIN POLARIZED INJECTION AND DISCUSSION 

Similarly to Ref. [Q] let us consider the case of a wide wave packet (approximated by a plane wave) initially injected into the 
lowest subband with spin up and an initial (group) velocity v g = yHko Y /m* along the y axis. In this case we find from the above 
equations 

(a~ l (t))= l-^y^sin 2 (Qf), (10) 

<u£(t)> = <o^(0> = 0. (11) 

To better understand the above result we should recall that, as shown in Ref. [[]], the corresponding expectation values of xn(t) 
and yH(t) are 

(xhO)) = ^sin 2 (Qf), (12) 
(yH(t)) = + ^-^- T [sin(2nf) -2Qf]. (13) 

This oscillatory zitterbewegung shown in the above expectation values is proportional to the initial group wave length, a fact 
which offers way more favorable perspectives for the experimental detection of this effect than for free or weakly bound electrons. 
In the latter case the amplitude of the zitterbewegung is of order the free-electron Compton wave length and therefore by orders 
of magnitude smaller. For further details we refer to the discussion in Ref. [[]]. Note also that differently from the zitterbewegung 
following from the ordinary Rashba s-o coupling in a single band [Q|, for which the trembling motion is perpendicular to the 
direction of propagation, here we find an oscillatory contribution in (yH(t)) along the initial group velocity. From Eqs. < JT~o| ) 
and ( |T2| ) we can see that (cr^(f)) = 1 - 2ko y (xH(t)}. That is, the dynamics in spin and real spaces are coupled. To make this 
point even more explicit, let us calculate the time dependent expectation of the z component of the orbital angular momentum: 
l z H (t) = U'U Z <g> 1U = U^{xpy - yp x ) ® 1U. A straightforward calculation yields 

<&(*)> = j Sin2(Qf) = ^OyiXHit)) (14) 

Since s = focr/2 defines the spin angular momentum, we can immediately see that 

Mm + (i zH (t)) = n/2. (15) 

The above result is easily undertood when we recall that at t — we injected a spin-up wave with zero angular momentum (i.e., 
at t — j z — s z = h/2). The conservation law in ( p"5| ) is general (i.e., valid not just for the expectation values) and follows from 
the fact that [L + s z , < K] = 0. The constraint / ; = -s z together with the s-o coupling strongly tie the dynamics of the system thus 
leading to the zitterbewegung in both the spin and real variables. 
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6. SUMMARY 



In summary, we have derived expressions for the zitterbewegung in the spin variables due to the inter-subband-induced s-o 
coupling present in symmetric wells with two subbands. For the simple case of an injected spin-up plane wave, we explicitly 
show that the zitterbewegung in both the spin and real spaces are intrinsically linked. This follows from both the s-o coupling 
and the conservation of the z component of the total angular momentum in our system. 
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